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We present a Green's function approach based on a LCAO 
scheme to compute the elastic propagation of electrons in- 
jected from a STM tip into a metallic film. The obtained 
2D current distributions in real and reciprocal space furnish 
a good representation of the elastic component of Ballistic 
Electron Emission Microscopy (BEEM) currents. Since this 
component accurately approximates the total current in the 
near threshold region, this procedure allows - in contrast to 
prior analyses - to take into account effects of the metal band 
structure in the modeling of these experiments. The Au band 
structure, and in particular its gaps appearing in the [111] and 
[100] directions, provides a good explanation for the previ- 
ously irreconcilable results of nanometric resolution and simi- 
larity of BEEM spectra on both Au/Si(lll) and Au/Si(100). 
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I. INTRODUCTION 

Ballistic Electron Emission Microscopy (BEEM)0i is 
a new technique based on the Scanning Tunneling Mi- 
croscope (STM). It has been primarily designed for the 
study of buried metal-semiconductor interfaces, in par- 
ticular for the investigation of the Schottky barrier. The 
experimental setup consists of a STM injecting current 
in a metallic film deposited on a semiconductor material. 
After propagation through the metal, a fraction of these 
electrons still has sufficient energy to surpass the Schot- 
tky barrier and may enter into the semiconductor to be 
finally detected as BEEM current. Using the tunneling 
tip as a localized electron source gives BEEM its unpar- 
alleled power to provide spatially resolved information on 
the buried interface, that can additionally be related to 
the surface topography via the simultaneously recorded 
tunneling current. 

The energy of the electrons contributing to the final 
BEEM current depends on the bias voltage between tip 
and metal, and is typically 1 to 10 eV above the Fermi 



energy in the metal. For energies close to the threshold 
voltage (the ones of primary interest for exclusive Schot- 
tky Barrier Height (SBH) investigations), thin metallic 
layers, and low temperature, the main contribution to the 
BEEM current stems from the elastic component, which 
are electrons not having suffered losses from electron- 
electron and/or electron-phonon interaction. Only in this 
limit can the technique be properly called ballisticB, and 
we shall concentrate in this paper on the propagation of 
such electrons and their contribution to the BEEM cur- 
rent. 

The Au/Si interface has been one of the first systems 
investigated by BEEM and has also proven to be a con- 
troversial one. Two important and apparently irrecon- 
cilable results have been (i) the elegant demonstration of 
nanometric resolution (« 15A) after propagation through 
films of as much as 100 — 150Aa, and (ii) the very sim- 
ilar results. obtained for the interfaces with Si(lll) and 
Si(100)Eltffl despite the strongly different k-space distri- 
bution of the projected Conduction Band Minima (CBM) 
available for injection of electrons into the semiconduc- 
tor. 

The theoretical analysis of BEEM data is usually based 
on the so-called four step modela: 1) tunneling from 
the tip to the surface, 2) propagation of the electrons 
through the metallic layer, 3) injection into the semi- 
conductor and 4) effects associated to various current- 
changing processes in the semiconductor (e.g. impact 
ionization and/or creation of secondary electrons). The 
last process becomes, however, only important at rather 
high energies and can be safely neglected when concen- 
trating on the near threshold regime. The standard 
model used to extract SBH's from the experiment is a 
simple quasi ID semiclassical approximation using pla- 
nar tunneling theory and WKB approximation for step 
1), free electron propagation plus simple exponential at- 
tenuation for step 2) and the QM transmission coefficient 
derived for a ID step potential for step 3). This predicts 
a 5/2-power law for the onset of the BEEM current above 
the threshold energy and has been used extensively for 
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the fitting of experimental dataH. 

Despite its obvious merits and especially its highly wel- 
comed simplicity, the standard model fails for more than 
a mere qualitative explanation of experimental data. In 
particular, the above mentioned results for the Au/Si sys- 
tem are not comprehensible within the framework of this 
crude approximation. As has been pointed out recentlyQ, 
the major error leading to these discrepancies for Au/Si 
is the complete neglect of band structure effects inside 
the metal film. The free electron treatment of the ex- 
isting model predicts an electron beam propagating in 
normal direction through the metal fibna., grown in [111] 
direction on either Si(lll) or Si(100)tffl. On the other- 
hand , it is a wellJaiown that Au shows a band gap in 
just this directiont2l, i.e. the Au band structure clearly 
forbids electronic propagation normally through the film. 

Unfortunately, the inclusion of band structure ef- 
fects into a theoretical model requires a much higher 
level of sophistication than for the previous simple ap- 
proximation based on E-space Monte-Carlo and various 
parametrized processes. We therefore present in this pa- 
per a Linear Combination of Atomic Orbitals (LCAO) 
scheme for the fully quantum mechanical computation 
of the elastic component of the BEEM current: step 
1) is included by coupling the tip to a semi- infinite 
crystal via the Keldysh Green's function formalismEJ, 
whereas for step 2) we can take the Au band structure 
fully into account via the corresponding Slater-Koster 
parameteral3'tj. This concept allows to compute the 
electron current distributions in real and reciprocal space 
in any layer inside the metal film. 

We will show in this paper, that already the detailed 
analysis of these 2D distributions immediately resolves 
the previously puzzling results for Au/Si as a pure con- 
sequence of the band gap enforced sideward propagation 
of the electron beam inside the metal film. This side- 
ward beam is highly focused thus explaining the obtained 
nanometric resolution, but it is also dominated by elec- 
trons with rather high k|| -momentum parallel to the in- 
terface, which allows similar injection into the projected 
Si CBM for both Si(lll) and Si(100) interfaces. Even 
though our proposed model is a little more demanding 
than the old simple one both from a theoretical and a 
computational point of view, we believe that it allows 
to gain deeper insight into the physics involved and that 
it represents a step towards a more careful treatment of 
the BEEM process which ultimately will enable a much 
better use of the obtained experimental data. 

II. THEORY 

Our model is based on a Hamiltonian written in a 
LCAO basis: 

H = Ht+Hs + Hj, (2.1) 
where Ht = e a n a + Tupc^cp defines the tip, 



Hs = J2 e i™i + J2Tijc\cj designates the metal sample 
and Hi = T am c' a c m describes the coupling between 
the tip and the surface in terms of a hopping matrix, 
T am (n a , cjj, and c a are number, creation and destruc- 
tion operators defined in the usual way). Note, that 
greek indices indicate tip sites, whereas latin indices cor- 
respond to sites inside the sample. The tight-binding pa- 
rameters Tij allow us to fully take into account the sam- 
ple band structure and can be obtained within the em- 
pirical tight-binding framework by a fitting to ab-initio 
band structuresEl Tabulated values exist for all elemen- 
tal crystals, in particular for Au as used in the present 
papercJ. The hopping matrix between tip and sample 
T am is also expressed as a function of the different atomic 
orbitals of atom a in the tip and atom j in the sample. 
Since BEEM experiments operate at high bias and rather 
low current, where the tip-surface distance is already of 
the order of 5 — 10 A, only tunneling between s-orbitals 
needs to be considered without loss of generality. The 
actual matrix element for this tunneling between the s- 
orbitals is simply described by an exponential function, 
as derived from the WKB approximations. It is impor- 
tant that T am is localized in a small region of m atoms 
close to the tip. 

A Green's function formalism presents the important 
advantage of being free of any adjustable parameters in 
the strictly elastic limit, where only an arbitrarily small 
positive imaginary part rj is added to the energy, E, nec- 
essary to ensure mathematical convergence: 

G R (E) = £. (2.2) 
E - H + it] 

Although in this paper we are not interested in discussing 
inelastic effects, which will form part of a forthcoming pa- 
per, it is worth mentioning that a finite r\ can be used 
to give some attenuation to the wavefield, mimicking in- 
elastic effects that draw away current from the BEEM 
experiment. We shall see in the discussion of the k-space 
current distributions that a finite r\ effectively defines a 
coherence region of size approximately given by A c w . 
Inside this region, quantum mechanical effects are im- 
portant, but outside of it quantum coherence is lost and 
intensities rather than amplitudes should be added to 
compute the final wavefield. 

The system under investigation is out of equilibrium 
as soon as a bias between tip and sample is appliedE-3. 
In order to retain the ab-initio advantage of a Green's 
function formalism, but still to be able to couple the 
tip to the-sample, it is convenient to use the Keldysh 
technique^, which essentially represents the generaliza- 
tion of many-bodsz-. Green's function theory to systems 
out of equilibriumtj. Within this formalism, the current 
between two sites i and j in the sample can be written 
as: 

Ja =^J Tr&jiG+r - G+p}dE. (2.3) 
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The matrices G -1 are non-equilibrium Keldysh Green's 
functions that can be calculated in terms of the retarded 
and advanced Green's functions G R /G A of the whole in- 
teracting tip-sample system. We are interested only in 
the elastic component of the current, 



G 



(I + G Yi )g (/ + EG), 



(2.4) 



where T, R ' A is the coupling between tip and sample that 
can be build up from the hopping matrices, T a j, link- 
ing atoms in both subsystems, while g^ refers to the 
Keldysh Green's functions of the uncoupled parts before 
any interaction has been switched on. 

The retarded and advanced Green's functions of the 
interacting system can be further obtained from a Dyson- 
like equation that uses the Green's functions of the un- 
coupled parts of the system, g R /g A , and the coupling 
between the surface and the tip, E fl,A : 



gR,A = 



R,A 



~R,A-^R,AqR,A 



(2.5) 



After some algebra, it is shown in Appendix A that 
the current between two sites i and j in the metal can be 
obtained, in the lowest order of perturbation theory with 
respect to the coupling between tip and sample, from the 
following formula: 



Jij(V) 
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dE, (2.6) 



where the integration is performed between the Schottky 
barrier (eV a ) and the voltage (eV) applied between the 
tip and the sample. We notice that this equation shows 
spectroscopic sensitivity because the Green's functions, 



gfm , an d the hopping matrices, T ma , depend on the en- 
ergy, and also through the upper limit in the integral. 
The summation runs over tunneling active atoms in the 
tip (a, (3) and the sample (m, n), p a p is the density of 
states matrix: (g A p — g R g) — 2i{ip a p. The trace denotes 
summation over the atomic orbitals forming the LCAO 
basis. 

Note, that the detour via the Keldysh technique has 
enabled us to arrive at an expression that only in- 
cludes the (tabulated) tight-binding parameters in- 
side the sample, the hopping elements T ma /Ts n between 
tip and sample, and the retarded and advanced equi- 
librium Green's functions g R /g A of the isolated tip and 
sample. 

We shall describe the injection of electrons from the 
metal into the semiconductor simply by matching the 
corresponding states. Therefore, we shall need to com- 
pute the current distribution at the interface in recip- 
rocal space. Working in a similar way as above for the 
real space, we find the following expression for the cur- 
rent distribution in reciprocal space defined now between 
planes i and j: 



J«(k||,V) = ^ Tr [£i(k||)$f ro (k 



f mQ (k||)p^(k||)7> n (k||)ffA(k||) dE (2.7) 



where in this case the sum runs over layers in the tip (a, 
0) and sample (m, n) that are involved in the tunneling 
process. All quantities in eq. ( |2.7[ ) are the fc|| -Fourier 
transforms of the corresponding objects in formula ([2.6]). 
The total current between planes i and j can be obtained 
by summing the current distribution Jy(k|i,V) as ob- 
tained in eq. (2/7) over the 2D Brillouin zone. 

It is worth mentioning, that even though the total re- 
sulting current at the given STM bias V is obtained as 
an integral in energies down to the SBH, the dominant 
contribution will be just the one at V itself, due to the 
strong exponential dependence of the tunneling coupling 
terms T a j with energy. This shall allow us to present 
qualitative results displaying only this main contribution 
(calculated at the highest energy, eV) , although we have 
additionally verified that essentially the same effects oc- 
cur at lower energi es. 

^From equation (2J3) it is clear that the main objects 



of interest that still need to be calculated are the retarded 
and advanced equilibrium Green's functions g R m / ' g A ti for 
a metallic surface. We_compute these quantities using 
a decimation techniqueOcJ, which involves an iterative 
process allowing in each step to double in size an existing 
slab (the so-called "super layers" in Figure |l|a). Assum- 
ing perfect periodicity in both directions parallel to the 
surface one can easily calculate gj^ fl (k|| , E), the retarded 
Green's function of the surface superlayer which contains 
the tunnel-active layers m, n. The layer doubling of the 
slab is then repeated until its size is so large, that the 
two surfaces are effectively decoupled. Hence, we obtain 
the Green's function projected at the surface, ffii (k|| , E), 
and the transfer matrix of the systemll3. With these two 
quantities it is possible to calculate the Green's function 
g^(k||,l?), representing the propagation of an electron 
of momentum k|| and energy E from the surface super- 
layer 1 up to the layer j inside the semi-infinite crystal. 
The corresponding Green's function in real space is sub- 
sequently obtained by Fourier transforming: 



kii 



$0*11,25) 



,iki|r,| 



(2.8) 



wherej— the summation is performed over N special 
pointsEil in the 2D Brillouin zone with respective weights 
tWfc.. . Finally, it is straightforward, to obtain the corre- 
sponding advanced Green's function by transposing and 
conjugating: 



g A (E*)=[gf 1 (E)]' 



(2.9) 



Equations (2.6)-(2.£) define our basic approach to the 
BEEM problem. 
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FIG. 1. Description of the geometry of the system, a) lat- 
eral view: tip coupled to the semi-infinite crystal, which is 
composed of superlayers used for the decimation; b) Au(lll) 
surface depicting the two tip positions treated in the text: Tl 
directly on top of one atom, T2 symmetrical bridge site be- 
tween two atoms (atoms in the surface to which tunneling is 
considered are whitened). 



leave the question of modeling the transmission factor for 
a forthcoming publication^! We feel that more physical 
insight is gained by visually comparing how the current 
distribution in the metal matches the available states in 
the semiconductor, than by presenting the final BEEM 
current, which results from a summation over all states 
with right matching conditions inside the 2D Brillouin 
zone. Therefore, in this work we simply draw in the 2D 
Brillouin zone the Si CBM ellipses to aid the eye identify- 
ing the relevant regions. We stress that the final summa- 
tion tends to smooth details and, as the physically mea- 
surable observable is only the kn -integrated I(V) curve, 
this might be one reason that helps to explain why the 
previous simple models achieve in many cases decent fits 
to the data, even though it is quite obvious that the un- 
derlying 2D current distribution is grossly inconsistent 
with the metal band structure. 



III. REAL SPACE RESULTS 



Let us draw our attention to the full quantum- 
mechanical calculations in real space based on equation 
( |2.6|) and the decimation technique. The total procedure 
outlined in this section enables the computation of 2D 
current distributions at any layer inside the semi- infinite 
crystal. Using eq. (2.6), and calculating Jij(V) for all 
atoms i inside a given layer, with current contributions 
from all their respective neighbours j in the layer above, 
we arrive at the current distribution in real space. This 
simulates the BEEM current impinging on the semicon- 
ductor after propagation through a metallic film with a 
thickness corresponding to the chosen layer i. Results of 
this type will help us to understand which spatial reso- 
lution can be expected in a BEEM experiment. Alterna- 
tively, the calculation of Jy (k|| , V) according to eq. (2.7) 
provides the current distribution between layers i and j 
inside the 2D projected Brillouin zone. Since ultimately 
the electrons need not only to have sufficient energy to 
enter the semiconductor as BEEM current, but also to 
have the corresponding k|| to get into the Si CBM, this 
type of plot will permit us to investigate the onset and 
strength of the BEEM current for a certain type of inter- 
face. 

When concentrating on the near threshold region, i.e. 
injection into the bottom of the Si conduction bands, 
the latter may be approximated by free electron-like 
paraboloids with appropriate effective massesE£l. In the 
2D projection on either the Au/Si(lll) or Au/Si(100) 
interface this will then lead to a number of ellipses in- 
side the Brillouin zone, representing areas through_which 
transmission into the semiconductor is possibleEj. To 
arrive at quantitative results, additionally a QM trans- 
mission factor for each of the matched states has subse- 
quently to be applied. However, in the present paper we 
shall be concerned rather with qualitative results appar- 
ent from the 2D distributions themselves, and may thus 



We apply our quantum-mechanical formalism to de- 
scribe the elastic propagation of electrons through a 
[lll]-oriented Au metallic layer. Due to the large lattice 
mismatch between Au (a,f C c = 4.08A) and Si (ddiamond = 
5.43A)E3, the formed interface is of rather poor quality. 
Although the growth of Au layers on Si has not yet been 
studied extensively, there is evidence from Low Energy 
Electron Diffraction and Auger studies, that suggest that 
Au films growing either on Si(lll)El or on Si(100)EI form 
crystals oriented in the [111] direction after the first four 
or five layers, which typically display heavy disorder, are 
completed. From the nature of our results, that show 
the gradual build-up of Bloch waves as a function of the 
distance to the surface, we deduce that the likely dis- 
order contained in four or five layers immediate to the 
interface, does not provoke the automatic destruction of 
the formed pattern. We therefore take the calculated 2D 
distributions in any given layer of the ideal semi-infinite 
crystal as a good approximation to the actual current ar- 
riving at the Si semiconductor after propagation through 
a Au film of corresponding thickness and passage through 
the disordered interface region. 

Before discussing our full quantum-mechanical results, 
it is worth analyzing the semiclassical limit, which is 
based on Kostcr'«-|approximation to compute the bulk 
Green's functions^ and an appropriate symmetrization 
of waves reflected at the surfaced. The central point 
in this semiclassical approach mainly valid for thick lay- 
ers is that it admits a simple geometrical interpretation 
for the bulk propagators, relating them to the shape of 
the constant energy surface (Fermi surface) of the metal 
band structure: applying a stationary phase condition all 
k-values contributing to the final Green's function may 
be neglected except one particular k c3: 
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surface, where the principal curvatures 




FIG. 2. Shape of the Au constant energy surface at leV 
above Ef- The shaded ellipses mark the regions contributing 
most to the projected 2D current distributions. The inset 
shows a schematic description for the relation between a real 
space propagation direction (Rj — Ri), and the relevant wave 
vector ko defined in the semiclassical approximation. 



R {E) = 



E 

k 



E - E(k) 



.ko(Rj-H-l) 



>kot?>ko(i) 



(3.1) 



dk'f dk'2 



where 



0^2, J-^f are the two principal curvatures associ- 
ated with the constant energy surface and the eigenvec- 
tors (j), </>k (l) are obtained diagonalizing the bulk 
tight-binding Hamiltonian H. 

For a given propagation direction (Rj — Ri), ko is 
the vector linking the k-space origin to the point on the 
constant energy surface (defined by the energy E) whose 
perpendicular to the surface is parallel to the initially 
chosen direction in real space (see inset in Fig. ||). This 
vector perpendicular to the constant energy surface at 
ko represents the group velocity of electrons being prop- 
agated in the metal, v g (ko) = VkE(k), that needs not 
to be parallel to ko when the metal band structure devi- 
ates from the free electron limit. Eq. (3.1) reveals that 



the current intensity is strongly enhanced in directions 
corresponding to rather flat parts of the constant energy 
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the denominator tend towards zero (note, that at the ex- 
tremal points, where the second derivatives are exactly 
zero, the next term in Koster's expansion is needed). 
Hence, already a simple visual analysis of the constant 
energy surface's shape may provide an intuitive under- 
standing where major ballistic current contributions are 
to be expected when including band structure effects. 

As pointed out previously^, for the analysis of electron 
propagation in Au layers at typical BEEM energies, it 
is crucial to notice that for voltages larger than w 1.0V 
necks develop in the [100] directions, which are similar 
to the well known necks in the [111] directions, cf. Fig. 
^. Identifying the flat parts on this surface and their 
corresponding real space propagation directions, we were 
able to make a qualitative prediction of the current dis- 
tributions in real and reciprocal space after propagation 
through a thick Au filmEl. 

In Fig. |2|, we have drawn the Au constant energy sur- 
face for 1 eV above Ep , and the regions (dashed areas) 
giving the major contribution to the semiclassical BEEM 
current in k-space. One of these areas is also shown in 
the figure inset: the reason why the k-space current dis- 
tribution maximum is associated with the A-region, and 
not with the apparently equivalent A'-region, is their dif- 
ferent orientation w.r.t. the (11 Indirection (see. Figs. 6 
and 8). It should be kept in mind that the dashed areas 
of Fig. define the momentum and velocity components 
of the electrons contributing most to the BEEM-current: 
the group velocity defining the semiclassical propagation 
in real space. Since most BEEM modeling up to now has 
been performed using E-space Monte-Carlo techniques, 
it may be interesting to notice, that apart from the neck 
regions where no propagation is possible, the Au constant 
energy surface remains spherical to a good approximation 
(Fig. g): the group velocity v g (k ) perpendicular to the 
surface and the connecting k- vector k do not diverge by 
more than « 20 — 30° for most directions of interest out- 
side the necks. When properly suppressing the forbidden 
propagation directions corresponding to the necks, a sim- 
ple application of k-space Monte-Carlo simulation should 
thus be possible, resulting in analogous results to ours in 
the purely ballistic limit. This allows the simulation of 
inelastic interactions in an approximate but simple way 
through the inclusion of the appropriate k-space scat- 
tering cross sections, reaching the interesting conclusion 
that those effects do not significatpely degrade the elastic 
predicted resolution in real spacec-i 

The presented semiclassical approach helps to better 
understand the numerical results in terms of the geomet- 
rical interpretation. However, we have also at hand the 
quantum mechanical decimation technique for the com- 
putation of the Green's functions, which allows us to pre- 
cisely calculate these distributions and also not only for 
thick films, but for any thickness. Since we are only in- 
terested in the qualitative evolution of injected current 
in a [lll]-oriented Au film, we start with a simplified 
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tunneling geometry, in which hopping is only permitted 
between one tip atom and one sample atom 1 position- 
ing the tip directly on top of the latter at 5.0 A distance, 
cf. Fig. [j]b. More realistic tunnel geometries allowing 
hopping to more atoms lead to essentially the same re- 
sults: one advantage of properly taking the metal band 
structure into account is that the results are not that 
much dependent on the initial tunneling conditions as 
long as metallic films of a certain thickness are involved. 
This is a strong difference between our theory and the 
standard ballistic E-space Monte-Carlo approach. The 
k-space distribution arriving at the semiconductor under 
the assumption of free electron propagation is inappro- 
priate for both very thick and very thin layers: (i) it is 
physically unreasonable to assume that what impinges 
on the interface (after a long propagation through the 
metal) is just the initial tunneling distribution, and (ii) 
on the other hand, when the film is thin enough for the 
details of tunneling to become important, the crude pla- 
nar theory based on WKB should not be applied. On 
the contrary, in all the examples presented in this pa- 
per (where we have deliberately avoided ultra thin lay- 
ers), tunneling simply provides a starting configuration, 
whereas the final current distribution at the interface is 
dominated by the preferred propagation directions dic- 
tated by the metal band structure. In any case, should 
the initial tunneling distribution become important for 
some particular application, our formalism would permit 
a more realistic point of view for that part of the problem. 

This discussion is illustrated in Fig. |[ where the re- 
sulting real space current distribution after electron prop- 
agation through more and more layers is depicted. As 
electron energy we chose 1 eV above the Fermi level of 
the metal. 2971 special k-points were used inside the 
2D Au Brillouin zone to suppress possible aliasing effects 
in the involved Fourier transforms to real space. In the 
second layer, very near to the surface (Fig. |3|a), the cur- 
rent is still strongly concentrated in the atoms closest to 
the location in the surface layer in which the current was 
injected. This behaviour, where the current still prop- 
agates in all directions and where the resulting distri- 
bution is basically a consequence of the crystal geome- 
try and nearest neighbour hopping, can be found down 
to about the fifth layer. Then, however, an interesting 
change occurs, which is related to the gradual formation 
of a Bloch wave inside the crystal: propagation becomes 
only possible in directions in accordance with the metal 
band structure. The immediate consequence of the al- 
ready mentioned band gap of Au in the [111] direction is 
that the evolving beam has to make way sidewards. This 
becomes obvious in Fig. [pD, where the opening up of the 
injected beam can already be observed with very little 
current remaining in the central plane atoms lying in the 
[111] forward direction. 

A further effect of the band structure is the forma- 
tion of narrowly focused, Kossel-likc lines correspond- 
ing to preferred propagation directions. After propaga- 
tion through a larger number of layers, a steady state 
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FIG. 3. Real space BEEM current distributions for 
Au(lll) after injection of current in one atom in the first 
layer (located in the center of the drawn layer). The tip is 
positioned at 5.0A height on top of the active atom. Param- 
eters: V = Ep + leV, r\ = O.leV. Distribution in a) 2nd layer 
(2.35A), b) 10th layer (21.19A), and c) 25th layer (56.5lA). 
Each dot represents one atom in the corresponding layer and 
the grayscale indicates the amount of current passing through 
the atom: black for maximum current to white for zero cur- 
rent. 
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FIG. 4. Intensity profile along two perpendicular cuts 
through one of the Kossel-like lines (the locations of the cuts 
are marked in Fig. ^c by white lines). 

is reached and distributions for consecutive layers differ 
only by the natural spreading of the resulting triangle 
as the electrons follow these preferred propagation direc- 
tions of the band structure, cf. Fig. ||c. These directions 
have a polar angle of about 30° with respect to the normal 
direction, and are in -good agreement with our previous 
semiclassical analysis^ (see Fig. |^). Note, how different 
the resulting pattern in each layer is from the prediction 
of the simple, free electron model: there the propaga- 
tion would be solely dictated by the k|| -spreading of the 
initial tunneling distribution resulting in a current cone 
centered around the [111] direction. Consequently, the 
real space distribution in any layer parallel to the sur- 
face would be a filled circle (or rather a filled triangle, 
considering the symmetry of a fee (111) crystalEa) with 
the major current contributions directly in the middle 
corresponding to normal propagation through the metal 
film. 

The deeper the chosen layer, i.e. the larger the thick- 
ness of the experimental film, the more spread over a 
larger triangle the current would be. The simple applica- 
tion of the uncertainty principle to the tusaeling process 
plus a free electron model for the metalala predicted a 
BEEM resolution for relatively thick 100 - 150A Au films 
of at best w J450A. On the contrary, the experiment by 
Millikcn et al.13 finds typical resolutions for such films on 
both Si(lll) or Si(100) of about 15A. In that experiment, 
sharp Si02 steps are created on both Si substrate types 
and subsequently covered by w 150A Au films. BEEM 
images of the step riser, when the tip crosses from above 
a part corresponding to the impenetrable Si02 to a part, 
where BEEM injection into the Si is possible, give a di- 
rect reflection of the size of the electron beam at the in- 
terface. The typical 15 A found are absolutely incompat- 
ible with the standard ballistic, free particle propagation 
used in the simple model and are difficult to justify even 
if parametrized electron-electrpn/electron-phonon inter- 
action is taken into accountcm. On the other hand, we 



believe that the observed formation of narrowly focused 
Kossel-like lines as caused by the Au band structure may 
already explain the experimentally obtained nanometric 
resolution in the purely elastic limit. As can be seen 
in Fig. ||c, the width of these lines carrying the BEEM 
current is typically 3-4 atomic distances. This is better 
appreciated in Fig. ^, where cuts perpendicular through 
these lines at two different positions were performed and 
the intensity vs. location is depicted. The derived width 
of rs 10A is in very good agreement with the experi- 
mentally observed value and is clearly distinct from the 
value derived from free electron theory, which deviates 
by a complete order of magnitude. 

A naive interpretation of this result seems to imply 
that multiple images of interface objects should be ob- 
served superimposed, since several such Kossel-lines are 
present in our distribution. A tip scan would sweep each 
of these lines across the interface object, each time lead- 
ing to a focused signal. It should, however, be born 
in mind, that so far we have used unrealistically sym- 
metric conditions for the current injection process: the 
on top tip site used, allows the current pattern ta dis- 
play the full threefold symmetry of the fee crystaBl and 
hence predicts the existence of three equivalent Kossel- 
like lines. During the scan, less symmetric tip sites are 
on the contrary involved, reducing logically the resulting 
overall symmetry of the current distribution. To verify 
this, we repeated the calculation, this time using a bridge 
site for the tip and allowing tunneling to both closest 
atoms on the surface, cf. Fig. |]b. 

Several conclusions can be drawn from the result shown 
in Fig. |[ which is for the limiting case after propagation 
through already a thick film. First, as expected, the over- 
all symmetry of the pattern has reduced from threefold to 
twofold, with one of the three Kossel-like lines decreasing 
in intensity. Second and most importantly, the shape of 
the distribution has essentially not changed. This is the 
confirmation of the introductory statement made at the 
beginning of this section: the inclusion of band structure 
effects reduces considerably the influence of the initial 
tunneling distribution. The preferred directions of prop- 
agation are exclusively dictated by the band structure, 
whereas it is only the population of these directions with 
electrons that depends on the actual tunneling process. 
In this respect, we again observe in Fig. |s| the w 30° 
off-normal, sideward propagation along focused lines, of 
which now only the two corresponding to the chosen tun- 
nel symmetry are mainly populated. 

As we are interested in highlighting the relevant phys- 
ical phenomena, we have introduced convenient approx- 
imations: one atom tip, spherical s-wave tunneling and 
a perfect substrate without any tilt, nor defects on the 
surface. In turn, the model neglects a number of fac- 
tors that necessarily would result in a lower symmetry 
than the one existing in a realistic experimental setup. 
However, it is plausible to accept that different factors 
concur in the experiment to degrade symmetry, and will 
eventually lead to the selection of only one of the pos- 
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FIG. 5. Real space BEEM current distribution as in Fig. 
^3, but injecting current under lower symmetry conditions: 
the tip is located in a bridge site at a height of 5.C)A. 



siblc Kossel-like lines, whose dominant contribution will 
then be responsible for the observed BEEM image with 
its nanometric resolution. On the other hand, none of 
the experimental influences can change the metal band 
structure itself (as long as there are at least crystal grains 
of the order of 100A), so that the formation of focused 
Kossel-like lines will be very similar to the one seen in 
our idealized model calculations. 



IV. RECIPROCAL SPACE RESULTS 

In order to understand how the current distribution in 
the metal matches the available states in the semicon- 
ductor, we now proceed to calculate the 2D distributions 
in reciprocal space using expression (2.7). Fig. [j]a shows 
the result obtained inside the Au Brillouin zone, again 
for the symmetrical on top tip site, cf. Fig. ^Jd, and for 
an electron energy of E = Ep + leV, which is still near 
the threshold region for Au/Si. As a direct reflection of 
the Au [111] band gap, which projects onto F in the cen- 
ter of the hexagon, the distribution displays a ring-like 
structure, since there are no states to carry the current in 
the r vicinity. It is again interesting to notice, how differ- 
ent this pattern is from the free electron prediction, that 
would inject the main current with k|| = just into F. 
In reality, complete reflection occurs for tunnel electrons 
with k|| — 0, permitting only injection into the Au(lll) 
states of rather high k|| -momentum. 

We would like to stress that the use of the decimation 
technique, -ajhich was derived from renormalization group 
tcchniqueaij, permits the fully quantum mechanical cal- 
culation of both required Green's propagators g^ildti °f 
the semi-infinite slab. The resulting BEEM current is 
consequently also fully quantum mechanical. It is there- 
fore not surprising that the distribution shown in Fig. |^a 
possesses a sixfold symmetry, which relates to the sym- 
metry of the [111] projected density of states to which 





FIG. 6. Reciprocal space BEEM current distributions in 
the 2D interface Brillouin zone, a) sixfold symmetry inside 
the coherence region, r\ — 0.001 eV b) threefold symmetry 
far outside the coherence region, r\ = 0.1 eV (2971 special 
k-points used in the 2D Brillouin zone, E = Ef + 1 eV; black 
dots represent high current, white zero current. A quadratic 
grayscale has been applied to emphasize the current changes 
due the symmetry crossover). 

+k and — k states contribute equallyS. 

This is at variance with a semiclassical distribution, 
where only k-vectors representing propagation towards 
the interface must be considered. These waves are clas- 
sically separated from those that propagate in the op- 
posite direction and the total vector field has thus the 
same symmetry as the constant energy surface, which in 
turn reflects the symmetry of the lattice, i.e. threefold 
for [111] directions in a fee materiaEj. Hence, our pre- 
vious semiclassical approaclvvielded threefold symmetric 
BEEM current distributions!!]. 

As has already been mentioned in section II, the use 
of a finite self-energy 77 in the construction of the Green's 
functions (2.2) puts some limits on the coherence length 
of our quantum mechanical computation. In the purely 
elasti c lim it, corresponding to a theoretical rj — > in 
eq. (2.2), the obtained distribution would not change 
at all after propagation through an infinite number of 
layers. In the ideal crystal, there are simply no processes 
enforcing k|| -changes of the elastic electrons after their 
injection. Numerically, we are however forced to choose 
a small, but finite r\ to ensure convergence. The result 
presented in Fig. ||a was hence obtained for r\ = 0.001 
eV, which is very small compared to the energy E = 1 
eV considered, and corresponds to a coherence length of 
several thousand Angstr0m. No significant changes of 
the distribution can be observed during the propagation 
through many layers. Since in this case we are always 
well inside the coherence region, the symmetry is sixfold, 
as it should be for a quantum mechanical calculation. 

On the other hand, we can considerably reduce the 
coherence length by using a much larger value for the 
optical potential, say r\ = O.leV. There would then ex- 
ist a regime outside the coherence region where interfer- 
ence terms do not play a significant role any more. This 
should allow us to recover the semiclassical results. Fig. 
Ob shows the obtained distribution for such a large r] and 
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Ay along cut 

FIG. 7. Intensity profile along the high symmetry line 
K — r — M of Fig. ^. The solid line represents current far 
outside the coherence region (cf. Fig. ^b), the dotted line 
current inside of it (cf. Fig. Ha). The different 77's lead to 
a different amount of damping in the currents and hence dif- 
ferent absolute scales. To enable better visual comparison, 
the currents have therefore been renormalized to the same 
maximum value. 

for the 30th layer in the crystal, i.e. well outside the co- 
herence region. It is most gratifying that not only the 
expected threefold symmetry appears, but that also the 
actual pattern matches perfectly with our semiclassical 
prediction, that was then only derived schematically, cf. 
Fig. 3 of ref. ( |^). Moreover, it is interesting to notice 
how our Green's function calculation permits to repro- 
duce the involved change in symmetry: even for a large 77, 
propagation in the first layers is still inside the small co- 
herence region and thus quantum mechanical effects are 
observed to produce the corresponding sixfold symmetry. 
During the propagation through more and more layers, 
the quantum coherence is gradually lost as is the sym- 
metry, which progressively approaches its threefold limit. 
This allows to study the crossover between the quantum 
and the semiclassical domain and provides another ex- 
ample of how a quantum system, under the influence of 
friction, becomes classical by a decoherence processEJ. 

Note, that the change in symmetry implies also a 
change of the detailed current distribution itself: the 
maxima of the semiclassical pattern appear on the lines 
r — M, whereas the quantum mechanical maxima lie 
along the directions T — K. This is more obvious in 
Fig. j^, where an intensity cut along the symmetry lines 
K — r — M is presented: the relative weight between cur- 
rents along K — r and T — M respectively is inversed 
between the two types of calculations, which is a re- 
flection of the symmetry change. It is interesting to 
notice that the evolution from the quantum symmetry 
to the classical one is gradually build up, and no sud- 
den change between both regimes is found. Our calcu- 
lations showed that the observed difference in reciprocal 
space does not significantly affect the beams in real space 
(where the symmetry must always be threefold), but it 
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FIG. 8. Matching of the obtained semiclassical k-space 
BEEM current distribution with the available CBM states 
in the semiconductor, which are approximated by projected 
parabolic bands. Current can only enter into the Si through 
the area enclosed by the ellipses [E — Ef + leV). a) 
Au/Si(lll), Esbh = 0.85eV b) Au/Si(100), E S bh = 0.82eV. 
Note, that the different lattice parameters of Au and Si re- 
quire remapping of the Si ellipses inside the larger Au Bril- 
louin zone. T—M corresponds to the [101] direction in k-space 

could nevertheless in principle affect the I(V) current in- 
jected through the projected ellipses into the semicon- 
ductor. However, because of the gradual crossover, no 
dramatic effects are to be expected, unless one could ex- 
perimentally break the time-reversal symmetry suddenly 
(e.g. by application of a magnetic field). 

As a final point we address the actual transmission of 
the current into the semiconductor states. The chosen 
energy E = Ep + leV is still in the near threshold region 
(E S BHjUu/Si{Ul)) « 0.85eV, E SBH (Au/Si(l00)) « 
0.82eVrj) justifying to approximate the Si CBM by nearly 
free electron ellipsoids. Since Si is an indirect semi- 
conductor, these ellipsoids are not located around T p, 
the Brillouin zone, but ca. 85% in direction T — Xt3 
and hence, [ttroject differently onto the [111] and [100] 
directionsnCJ. The resulting ellipses for both orienta- 
tions are drawn in Fig. ||a and ||b. Note, that since we 
are calculating inside the larger Au Brillouin zone, those 
ellipses belonging to higher Si Brillouin zones, but still 
within the first Au one, have to be considered for cur- 
rent injection as well. The principal difference for both 
Au/Si(lll) and Au/Si(100) interfaces is that in the lat- 
ter case Si ellipsoids project directly onto T, which is 
where the simple free electron picture puts the maxi- 
mum current. Since very little current carried by elec- 
trons with high ky -momentum is predicted in this type 
of model, a considerable difference in the onset and ab- 
solute magnitude of the BEEM curxcnt was originally 
anticipated between both interfaces™. On the contrary, 
the actual, experimentally observed spectra were highly 
similarErElu. Since this result was irreconcilable with the 
ballistic free electron theory, a variety of processes rang- 
ing from more isotropic tunnel distributions!] over strong 
elastic electron-electron inlfiractionB to ky -violation at 
the non-epitactic interfaces^ were proposed in an exten- 
sive debate in the literature. All these processes aimed 







at providing k|| -momentum to the current distribution 
to allow injection into the off-normal [111] ellipses, but 
ran into considerable trouble through the simultaneous, 
inherent loss of resolution!! 

On the other hand, the ring-like current pattern pre- 
dicted by our band structure calculation is dominated 
by high k|| -momenta, (cf. Fig. 8) but without loss of 
resolution, as we have seen in the last section. The 
total reflection of the injected electrons around F due 
to the projected Au band gap renders the existence or 
non-existence of a Si ellipse in that region completely ir- 
relevant, since no current can enter through it into the 
semiconductor in any case, and the current is therefore 
forced to enter via any of the off-normal ellipses for both 
Au/Si(lll) (Fig. |a) and Au/Si(100) (Fig. |b). The 
ring-like distribution may in zeroth order be well de- 
scribed as azimuthally symmetric, making it highly plau- 
sible why the BEEM current shows similar onset and 
magnitude: the exact location of the off-normal ellipses 
plays only a minor role. 

Finally, we quantify this result by comparing the actual 
current injected into the (100) and (111) silicon orienta- 
tions. As mentioned in the beginning, we shall not get 
at this stage into the separate problem of computing a 
suitable transmission coefficient to keep our discussion as 
simple as possible. Hence, we simply use T = 1, which is 
justified because we shall only present the ratio between 
the current injected into both silicon faces Jioo/^m, 
where factors not very sensitive to a particular silicon ori- 
entation cancel out. In particular, this is the case for the 
transmission coefficient T(E, k|| ) as can easily be_ checked 
with the crude ID step-barrier approximationo. Simi- 
larly, we do not need to consider inelastic effects, that 
are anticipated to affect equally electrons approaching 
the Si(100) or the Si(lll) surface, and should in general 
be small when considering thin metal films. The ratio 
shown in Fig. ^ compares well with the experimental 
ratio obtained from refs. ( |32| , |33| ), reproducing not only 
the correct order of magnitude, but also the overall trend 
with energy. The largest difference of «30% is seen in the 
low voltage region where experimental currents are weak 
and difficult to measure, and where the total current de- 
pends most sensitively on the exact value of the SBH for 
both orientations, whereas we have made no effort in fit- 
ting these values, but have simply taken averages over 
the big number of existing experimental values (which 
differ easily by 0.1 eV). It is hence not surprising that the 
agreement in the lowest 0.2 eV above the SBH is not as 
convincing as it is at the higher voltages. Taking finally 
into account the crudeness of the applied model where we 
have deliberately avoided the fitting of any parameter to 
focus on physical insight, we believe that the presented 
ring- like distribution provides a satisfactory and intuitive 
explanation of the observed effect. 
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FIG. 9. Ratio of current injected in Si(100) and Si(lll) af- 
ter propagation through thin Au films of approximately equal 
thickness (100 A Au/Si(100) and 75 A Au/Si(lll), chosen to 
compare with available experimental data). Elastic electron 
propagation and T=l are used in then-theory (dashed line); 
experimental data taken from refs. (Olj) (solid line). 

V. CONCLUSIONS 

In conclusion, we have presented a theoretical model 
that allows the calculation of the elastic contribution to 
the BEEM current in the near threshold range. Based on 
a LCAO scheme and a Keldysh Green's function formal- 
ism, 2D distributions in real and reciprocal space may be 
computed in any layer of a semi-infinite crystal after cur- 
rent is injected from a tip atom. Subsequent matching 
with the semiconductor states permits to extract the ac- 
tual BEEM current, although the main emphasis in the 
present paper has been on the qualitative understanding 
that can already be obtained from the 2D patterns them- 
selves. The model allows for the first time to fully take 
into account the influence of the metal band structure in 
the BEEM process, and possesses the further advantage 
of being free of any adjustable parameter in the strictly 
elastic limit. 

The application to the system Au/Si shows a variety 
of consequences of the Au band gap in the [111] direc- 
tion. In real space, the formation of narrowly focused 
Kossel-like lines and a sideward beam propagation is ob- 
served, which may explain the experimentally obtained 
nanometric resolution. In reciprocal space, a symmetry 
change between quantum mechanical and semi-classical 
regime can be related to the gradual breaking of quantum 
coherence. The actual pattern in k-space has a ring-like 
shape, which calls for current injection via the off-normal 
Si ellipses for both Au/Si(lll) and Au/Si(100). Hence, 
the sole inclusion of band structure effects achieves to 
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explain the nanometric resolution and the similarity of 
BEEM spectra on both Si orientations on the same foot- 
ing in the purely elastic limit - without any adjustable 
parameter and without the necessity to invoke any fur- 
ther scattering process. 



and 
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APPENDIX A 



The term (gf~ — gtj~) m (5.4) gives a zero contribu- 
tion to the current because it corresponds to the current 
between sites i and j inside the metal in the absence of 
tip-sample coupling, which obviously must be zero. 

Let us work now with the term associated to the tip, 
Ut- This term can be simplified by using the well known 
relation between advanced and retarded Green's func- 
tions: G A m = [G mn )^- The real matrix T ma is just equal 
to Hence: 



We start from formula (|2.3|) and the Keldysh Green's 



functions defined in (2.4). The coupling term S ' can 



be expressed as a function of the real hopping matrices, 
T Qm , that link tunneling active atoms in the tip (a) with 
the corresponding ones in the sample (m). 
Using equation 



G 



+- 

ji 
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am 



, Gf and Gf- can be written as: 



\_GfaTam9mi + 9jm TmaG ai 



/ j G j m T ma g a p Tp n G ni 



(5.1) 
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C H T — T 1 



(5.2) 



The Keldysh Green's function g H of the uncoupled 
(and hence equilibrium) system can be expressed as a 
function of the advanced and retarded Green's function 
of the uncoupled parts g R /g A and the Fermi distribution 
functions of the tip, fr(E), and sample, fs(E): 



(5.3) 
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Using these last equations we obtain: 

+" -J- f T U T + fsU s 



G t)~ 



9ji ~9 V 

where we have defined the auxiliary variables 



(5.4) 



5^ \g j m T ma 
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flR rp 
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(5.5) 
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(5.7) 



So Ut can be expressed as the real part of a matrix: 

U T — 2K 2J GfmTma(9£p — 9ap)T0nGni ( 5 - 8 ) 



ma/3n 



The retarded and advanced Green's functions for the 
interacting system can further be obtained from a Dyson- 
like equation that uses the Green's functions of the un- 
coupled parts of the system g R / g A and the coupling term 



E R / E A (equation |J): 
G R A = q RA - 



~R,Aj^R,AqR,A 



(5.9) 



If the values of the coupling matrix T am are much 
smaller than hopping terms inside the metal (as it is the 
case in tunneling conditions) we can work in the low- 
est order perturbation theory and approximate G R m and 

Gfi simply by g R m and g A it respectively. Moreover, we 
can further simplify the expression for Ut by relating 
(da/3 ~ 9 'a/s) w ith the density of states matrix at the tip 
p a p by the equation g A p - g R p = 2mp a p: 



U T = 4tt3 



^ ] \^ jmTmaPa/3Tp n g n i 



(5.10) 



ma/3n 



The term associated with the sample, Us-, can also be 
written as the re al pa rt of a matrix using similar argu- 
ments leading to ( |5.§| ): 



U S 



2^ y ] ^Gj a T a mgi m gj m r EmaG 



ma 



g jm^maG a i 



r R f n R 
u ja* amy mi 



(5.11) 



Using the Dyson-like equation for the Green's func- 
tions of the interacting system and again working in the 
lowest order of perturbation theory: 
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(5.12) 



(5.13) 



The term associated with the sample Us can then be 
expressed as: 



maf3n 

4tt9 ^ (gfmTmaPapTfln 
motftn 



(5.14) 



It is easy to demonstrate that the first term of this 
equation is zero because it is the real part of the difference 
of two magnitudes, one being the complex conjugate of 
the other. Hence Us can be shortened to: 



Us = -4tt3 ^2 9jLTmaP a pTf3 n g, 



(5.15) 



ma/3n 



So the rest between and Gj- is finally written: 



ma.(3n 



(5.16) 



Assuming zero temperature, we obtain the current be- 
tween two sample sites i and j in the sample as an inte- 
gral over a window of energies ranging from the Schot- 
tky Barrier Height (eV ) up to the applied voltage of the 
imaginary part of the trace of a product of matrices: 



4e, 



f eV r - 

I Tr \ T ij9rm T rnaPa/3T/3 n gni 
JeV n _ . o_ L 



dE 



mctfln 



(5.17) 



This is equation ( |2.6| ), that can be considered as our 
starting point. 
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